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We study the effects of laser-induced Rashba-like spin-orbit coupling along the Bardeen-Cooper- 
Schrieffer-Bose-Einstein condensate (BCS-BEC) crossover of a Feshbach resonance for a two-spin- 
component Fermi gas. We calculate the condensate fraction in three and two dimensions and 
find that this quantity characterize the crossover better than other quantities, like the chemical 
potential or the pairing gap. By considering both the singlet and the triplet pairings, we calculate 
the condensate fraction and show that a large enough spin-orbit interaction enhances the singlet 
condensate fraction in the BCS side while suppressing it on the BEC side. 
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I. INTRODUCTION 

Over the past several years, the predicted crossover 
[l|-[3[ from the Bardeen-Cooper-Schrieffer (BCS) state of 
weakly bound Fermi pairs to the Bose-Einstein conden- 
sate (BEC) of molecular dimers has been observed by 
several experimental groups 041] ■ In two experiments 
@, @ the condensate fraction of Cooper pairs , which 
is directly related to the off-diagonal- long-range order of 
the two-body density matrix of fermions [ll[ [TH, has 
been studied with two hyperfmc component Fermi va- 
pors of 6 Li atoms in the BCS-BEC crossover. The exper- 
imental data are in quite good agreement with mean-field 
theoretical predictions [lj, [l4| and Monte-Carlo simula- 
tions [HI at zero temperature, while at finite temperature 
beyond- mean- field corrections are needed [l6j . Recently 
the condensate fraction in the BCS-BEC crossover for 
a two-dimensional (2D) Fermi gas [13], and for a three- 
spin-component Fermi gas with SU(3) symmetry (l8| . has 
been theoretically investigated. Remarkably, last year 2D 
degenerate Fermi gases were experimentally realized for 
ultra-cold atoms in a highly anisotropic disk-shaped po- 
tential (l9j . 

Quite recently, artificial spin-orbit coupling has been 
obtained in neutral bosonic systems [20(, where the 
strength of the coupling can be controlled optically, and it 
has been suggested that the same techniques can be used 
with ultracold fermions [U [22[ . These results have stim- 
ulated the theoretical investigation of spin-orbit effects 
with Rashba [23]| and Dresselhaus 24 \ terms in the BCS- 
BEC crossover [25|-|30(. In particular, very recently and 
independently, several authors have analyzed the evolu- 
tion from BCS to BEC superfluidity in the presence of 
spin-orbit coupling for a 3D uniform Fermi gas [2oT - |28| . 
and in the 2D case by use of a perturbative approach 
f29l ]. Nevertheless, those papers did not considered the 
condensate fraction of Fermi pairs. 

In the present paper we calculate the chemical poten- 
tial, the pairing gap, and the condensate fraction along 
the BCS-BEC crossover both in 3D and 2D as a function 
of spin-orbit coupling. We show that the two contribu- 
tions - i.e., those related to the singlet and triplet pairings 



- to the condensate fraction, separately, characterize the 
crossover better than the other quantities. Remarkably, 
a large enough spin-orbit interaction enhances the singlet 
contribution to the condensate fraction in the BCS side 
while suppressing it on the BEC side. The triplet con- 
tribution to the condensate grows by increasing the spin- 
orbit coupling and is larger close to the crossover. On the 
contrary, the chemical potential and the pair function ex- 
hibit no peculiarities along the crossover. Moreover, we 
find that, when the Rashba velocity becomes of the order 
of the Fermi velocity, there is a value for the dimension- 
less interaction strength y — l/(kpa s ), where &f denotes 
the Fermi linear momentum and a s the inter-atomic s- 
wave scattering length, for which the singlet condensate 
fraction no longer depends on spin-orbit coupling. This 
nodal point can be promoted as the real point of the 
crossover. What is observed in three dimensions occurs 
also in two dimensions where the nodal point occurs when 
the binding energy cb is almost equal to the Fermi en- 
ergy cf, i.e., eb ~ £f- In 2D the condensate fraction 
approaches the value 1 only for extremely large values of 
the scaled binding energy cb/cf- 



II. THE MODEL 

Let us consider the following Hamiltonian 
H = Hq + Hi , 



(1) 



where Hq is the single particle Hamiltonian in the pres- 
ence of Rashba and Dresselhaus terms [23|, [24| , namely 
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where vd and vr are, respectively, the Rashba and Dres- 
selhaus velocities; a x and a y denote the Pauli matri- 
ces in the x and y directions; and ^(k) is the spinor 
■0(k) = (ipf(k), ipi(k)) T . Hi is the interaction term given 
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by 
Hi 



j- ^(k+q)^(-k)^(-k' + q)V>t(kO , (3) 

kk'q 



where g > 0, which corresponds to attractive interac- 
tion. After defining the order parameter describing the 
particle pairs, A = (g/V) X}k(^4-( — k)-0-|-(k)), where V is 
the volume, at the mean field level we can decouple the 
interaction, finding 



Hl = V l -^~Y, (AW-k)Vt(k) + AV{(k)V-|(-k) 

(4) 

Introducing the following multispinor ^(k) = 
(■0-|-(k), k), ^4,(k), k)) T , one can resort to 
standard path integral formulation at finite temperature 
obtaining, within a saddle point approximation and after 
integrating over the fermions [31], the thermodynamic 
potential 

• y M ku k 

where j3 — l/(kgT), where kg denotes the Boltzmann 
constant and T the absolute temperature, and G" 1 is a 
matrix on the basis of 4 r (k) which reads 
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with 7(k) = hvn(ky + ik x ) + hvr>(k y — ik x ) and £k = 
H 2 k 2 /2m — a. After summing over the Matsubara fre- 
quencies [31j the thermodynamic potential becomes 



« = ^ - ^ E E mi + e ~^ (k) ) +E& 



(7) 
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w here ^(k) = V (g k - | 7 (k)|) 2 + \Af, E 2 (k) = 
V '(a + l7(k)|) 2 + |A| 2 , E 3 = -E x and E t = -E 2 . From 
the thermodynamic formula N — — we obtain the 
equation for the number of particles 
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The gap equation is, instead, given by 
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and finally the condensate number [12|, |13j, |32j reads 

W o =JV + JVi, (10) 



where 

N = ]T|(^(k)^(-k))f 
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is the singlet contribution to the condensate, with total 
spin 0, whereas 



Ni = ^|^ t (k)^(-k)}| 

k 
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is the triplet one, with total spin 1. 
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FIG. 1: Singlet (no, upper curves) and triplet (m, lower 
curves) condensate fractions of the 3D Fermi gas as functions 
of the dimensionless interaction strength y — l/{kpa s ) for 
different values of Rashba velocity, (vr/vf) 2 = (solid line), 
0.5 (long-dashed line), 1 (short-dashed line), 2 (dotted line), 
and 4 (dashed-dotted line). 

We are interested in the low temperature regime where 
the condensate is quite large. Quantitatively we can re- 
strict our study to the zero temperature limit, or, in 
three dimensions, when at least 2k bT <C A, with A 
now supposed to be a real number. The zero temper- 
ature limit is also mandatory for the two dimensional 
case. In the equations above we have therefore simply 
tanh( / 3.E i (k)/2) -> 1. 



III. THREE DIMENSIONS 

Let us consider first the three dimensional case. Here- 
after we proceed in the same spirit of Ref. [33j], gener- 
alizing the calculation including the spin-orbit coupling. 
After rescaling the momenta 



v / 2mA 
k = — z — q 
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FIG. 2: Singlet condensate fraction of the 3D 
Fermi gas as a function of (vr/vf) 2 for y = 
-2,-1.5,-1,-0.5,0,0.18,0.5,1,1.5,2 (corresponding to 
the curves from below). The solid thicker line is for y = 0. 
The dashed line is when y ~ 0.2, below this value all the 
curves have a minimum, while above it they always decrease. 



and summing in the continuum (^ k — > f rf 3 k) we 
get, for the number of particles 



N _ (2mA) 
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.3/2 
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with dimensionless parameters denned as follows 
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x\ = 2m 
X2 = 2m 
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In the continuum limit, due to the choice of a contact 
potential, the gap equation ^ diverges in the ultraviolet. 
After regularization 0] the gap equation reads 
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where a s is the s-wave scattering length between fermions 
with different spin component. In this way we get 
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FIG. 3: 3D Fermi gas: Chemical potential /i (upper panel) 
and pair function A (lower panel), both in units of £f, as 
a function of the dimensionless interaction strength y — 
l/(ii?a s ) for different values of (vr/vf) 2 (th 



e same ones as 



in Fig. \T\ with same line types). 

Finally the condensate densities are given by 

iff (x ,x 1 ,x 2 ) , 



(2mA) 3 / 2 
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(20) 



lff(x ,xi 



r=±l 
2 



x + rJx\q 2 + x\q. 



r 2„2 



(21) 



We can write also the gap and the chemical poten- 
tial in terms of the Fermi energy tp = H 2 k F /(2m) = 
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h 2 /(2m)(3n 2 n) 2 / 3 , therefore 

^=4(|) a/3 4Wi,* 3 )- 2/3 . (22) 
^ = 4(^) 2/3 X oI 3 N d (xo,x 1 ,x 2 )- 2 / 3 . (23) 



Finally the spin-orbit velocities can be written in terms 
of the Fermi velocity 

^^=(f) 2/3 , 1)2 /^o,x 1 ,x 2 )-/3 (24) 

We are now in the position to express the two contribu- 
tions to the condensate fraction 
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2n a _ 1 lff g (x 0: xi,x 2 ) 
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the chemical potential, Eq. (|2"3"|) . and the gap, Eq. 
in terms of the scattering parameter y, Eq. (|18l) . This 
is guaranteed, at least heuristically, by the fact that 
for any point in the space of dimensionless parameters, 
(xo,Xi,X2), there are single values for y, 2n s /n, n/ep 
and A/ep. For x\ = x 2 = 0, namely without spin-orbit 
couplings, we indeed recover previous analytic results re- 
ported in Ref. 33]. 

The results shown here in the figures are obtained fix- 
ing x\ = x 2 , namely when only Rashba {vp> = 0) or only 
Dresselhaus (vr = 0) are present. The other special case 
with vr = vd is actually less interesting since in that 
case the singlet condensate fraction seems to be always 
suppressed. On the contrary, with only Rashba term (or 
only Dresselhaus) we observe (see Fig. [J) that, turning 
on the spin-orbit coupling, the singlet condensate frac- 
tion increases in the BCS regime, whereas it decreases in 
the BEC regime. In particular, for vr > vp the singlet 
condensate fraction at y s» 0.2, slightly above the unitar- 
ity, which is 2no/n rs 0.7, does no longer depends on the 
spin-orbit coupling; see Fig. [2] On the left (BCS side) 
of this point the singlet condensation is improved by the 
spin-orbital interaction, whereas on the right (BEC side) 
this condensation is suppressed. The triplet condensate 
fraction, instead, decreases in both the BCS and BEC 
limits, whereas it is sizable close to the crossover, ex- 
hibiting a non-monotonic behavior. On the contrary, the 
chemical potential \i (upper panel of Fig. [3]) is shifted to- 
ward negative values in both the regimes, while the pair 
function A (lower panel of Fig. [3]) is enhanced both in the 
BCS side and in the BEC one, although, in the latter, the 
enhancement is less pronounced. These last quantities, 
therefore, unlike the condensate fraction, exhibit no pe- 
culiarities across the crossover. 



IV. TWO DIMENSIONS 



In two dimensions the regularization of the gap 
tion differs, with a bound state always present [34j 



equa- 
With 



cb as the binding energy, we have, therefore, 



E 



i 



i i 



g 2(£ k + M ) + e B ' 



(26) 



which, after rescaling the momenta as in Eq. and 
integrating, leads to 



1 m ( 2k 2 
~g = irf ln U + 11 ' 
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where A is the ultraviolet momentum cut-off. On the 
other hand Eq. ^ holds, where now the sum is over 
momenta in two dimensions, and, therefore, we have 
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FIG. 4: Singlet (no, upper curves) and triplet (m, lower 
curves) condensate fractions of the 2D Fermi gas as functions 
of the binding energy eb, in units of the Fermi energy ep, for 
different values of Rashba velocity, (vr/vf)' 2 = (solid line), 
0.5 (long-dashed line), 1 (short-dashed line), 2 (dotted line), 
and 4 (dashed-dotted line). 

From this expression we derive the binding energy 

~ir = um m w — i — 7 ' (30) 

A A-voo exp[i s (xo, xi, X2)/7rj — 1 

which actually does not depend on the cutoff since I g has 
a logarithmic divergence in the ultraviolet which cancels 
exactly the factor A 2 . In the absence of spin-orbit, x\ = 
x 2 = 0, we recover, in fact, the known result cr/A = 



0.0 
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FIG. 5: Singlet condensate fraction of the 2D 
Fermi g function of (vr/vp) 2 for es/tF = 

0.01,0.06,0.2,0.5,1,1.5,2,2.5 (corresponding to the curves 
from below). The dashed line is when tB/ep ~ 1, below this 
value all the curves have a minimum, whereas, above it, they 
always decrease. 



As in the three dimensional case, the quantities we 
consider are the following 



2n s _ l rf? s {x ,xi ) x 2 ) 
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as functions of £b/ € f where now ep — h 2 Trn/m is the 
Fermi energy in two dimensions. The integrals Iff and 
iff are the same ones as in Eqs. (|13p and (|21|l but defined 
in two dimensions (d 3 q — > d 2 q). 

Experimentally, the realization of a 2D system can be 
obtained by a strong harmonic confinement in one direc- 
tion, i.e., lo z 3> u) x , uj y ; therefore one can link the tunable 
3D scattering length a s to the two-body binding energy 
in 2D. Introducing, conveniently, the confining length 
l z = \J h /muj z , one finds, in fact, m(es/ftwz) ~ £ z /a s 
(for more details see Refs. [H, H^|). 

Again we focus our attention to the case with only 
Rashba (or, equivalently, only Dresselhaus) term, i.e., 
x\ = x 2 . Also in this case, as in the three-dimensional 
one, the spin-orbit produces interesting effects in the con- 
densate fraction, showing a nodal point at cb ~ e_F, set 
in when vr ~ vp, see Fig. |4] and Fig. [5l in the neigh- 
borhood of which the slope of the curve decreases. As a 
result, the spin-orbit coupling promotes the singlet con- 
densation in the BCS side and suppresses it on the BEC 
side. We observe that, contrary to the 3D case, in 2D the 
condensate fraction approaches the value of 1 only for an 
extremely large interaction strength, i.e. for €b/^f 3> 1- 
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FIG. 6: 2D Fermi gas: Chemical potential /i (upper panel) 
and pair function A (lower panel), both in units of ef, as a 
function of the binding energy en (in units of ep) for different 
values of (vr/vp) 2 (the same ones as in Fig. [3] with same line 
types). 



Again, the chemical potential \x is pushed toward more 
negative values (see the upper panel of Fig. [6]). For small 
spin-orbit coupling, and eb — > 0, one recovers the known 
perturbative result, /i ~ ep — mv\/2. The pairing gap 
A (lower panel of Fig. \5§ is, instead, increased by the 
Rashba spin-orbit interaction in the whole crossover. 



V. CONCLUSIONS 

We have studied the evolution of BCS superconduc- 
tors to BEC superfluids in the presence of an artificial 
spin-orbit coupling of Rashba and/or Dresselhaus type in 
two and three dimensions. We have shown that, unlike 
the chemical potential and the pairing gap which exhibit 
no particular behaviors at the crossover, the condensate 
fraction is very peculiar. The condensation of singlet 
pairs, in fact, is promoted by Rashba coupling in the BCS 
regime whereas it is suppressed in the BEC regime. In 
the middle, both in three and in two dimensions and for 
large enough Rashba spin-orbit coupling, there is a nodal 
point where the curves of the singlet condensate fraction 
cross each other, and, for this reason, this can be con- 
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sidered the putative point of the crossover. On the other 
hand, the triplet contribution to the condensate fraction 
has not a monotonic behavior as a function of the scat- 
tering parameter, swelling up close to the crossover. The 
full condensate fraction increases with the spin-orbit in- 
teraction. Because in our calculations we have used the 
mean-field theory, it is important to stress that Monte 
Carlo simulations have shown that, at zero temperature, 
beyond-mean-field effects are negligible in the BCS side 
of the BCS-BEC crossover, whereas they become relevant 
in the deep BEC side [H, HH . In conclusion, we think 



that our results can be of interest for future experiments 
with artificial gauge potentials in degenerate gases made 
of alkali-metal atoms. 



Acknowledgments 

L.D. thanks the International School for Advanced 
Studies, SISSA, Trieste, for hospitality and facilities dur- 
ing the completion of this work. 



D.M. Eagles, Phys. Rev. 186, 456 (1969). 
A.J. Leggett, in Modern Trends in the Theory of Con- [17 
densed Matter, p. 13, edited by A. Pekalski and J. Przys- [18 
tawa (Springer, Berlin, 1980). [19 
P. Noziers, S. Schmitt-Rink, J. Low Temp. Phys. 59, 195 
(1985). [20 
M. Greiner, C.A. Regal, and D.S. Jin, Nature (London) 
426, 537 (2003). [21 
C.A. Regal, M. Greiner, and D.S. Jin, Phys. Rev. Lett. 
92, 040403 (2004). [22 
J. Kinast, S.L. Hemmer, M.E. Gehm, A. Turlapov, J.E. 
Thomas, Phys. Rev. Lett. 92, 150402 (2004). [23 
M.W. Zwierlein et al., Phys. Rev. Lett. 92, 120403 

(2004) ; M.W. Zwierlein, C.H. Schunck, C.A. Stan, [24 
S.M.F. Raupach, W. Ketterle, Phys. Rev. Lett. 94, [25 
180401 (2005). [26 
C. Chin et al, Science 305, 1128 (2004); M. Bartenstein [27 
et al, Phys. Rev. Lett. 92, 203201 (2004). [28 
Y. Inada, M. Horikoshi, S. Nakajima, M. Kuwata- [29 
Gonokami, M. Ueda, and T. Mukaiyama, Phys. Rev. [30 
Lett. 101, 180406 (2008). [31 
C.N. Yang, Rev. Mod. Phys. 34, 694 (1962). 

O. Penrose, Phil. Mag. 42, 1373 (1951); O. Penrose and [32 
L. Onsager, Phys. Rev. 104, 576 (1956). 
C.E. Campbell, in Condensed Matter Theories, vol. 12, [33 
131 (Nova Science, New York, 1997). 

L. Salasnich, N. Manini, and A. Parola, Phys. Rev. A 72, [34 
023621 (2005). 

G. Ortiz and J. Dukelsky, Phys. Rev. A 72, 043611 [35 

(2005) . 

G. E. Astrakharchik, J. Boronat, J. Casulleras, and S. [36 
Giorgini, Phys. Rev. Lett. 93, 200404 (2004). 
Y. Ohashi and A. Griffin, Phys. Rev. A 72, 063606 
(2005); N. Fukushima, Y. Ohashi, E. Taylor, and A. Grif- 



fin, Phys. Rev. A 75, 033609 (2007). 

L. Salasnich, Phys. Rev. A 76, 015601 (2007). 

L. Salasnich, Phys. Rev. A 83, 033630 (2011). 

K. Martiyanov, V. Makhalov, and A. Turlapov, Phys. 

Rev. Lett. 105, 030404 (2010). 

Y.J. Lin, K. Jimenez- Garcia, and LB. Spielman, Nature 
471, 83 (2011). 

J. Dalibard, F. Gerbier, G. Juzeliunas, P. Ohberg, e- 
preprint la^iv:1008.5378H (2010). 

M. Chapman and C.A.R. Sa de Melo, Nature 471, 41 
(2011). 

Y.A. Bychkov and E.I. Rashba, J. Phys. C 17, 6029 
(1984). 

G. Dresselhaus, Phys. Rev. 10 0, 580 (1955). 

M. Gong, S. Tewa ri, C. Zhang, larXiv:1105. 17961 
Z-Q. Yu, H. Zhai. larXiv:1105.2250l 

H. Hu, L. Jiang, X-J. Liu, H. P u,larXiv:1105.2488l 
Li Han and C.A.R. Sa de Melo, larXiv:1106.3613l 

G. Chen, M. Gong, and C. Z hang, | arXiv: 110 7.26271 
M. Burrello, A. Trombettoni, arXiv:1108.0839] 

H. T.C. Stoof, B.M. Dennis, and K. Gubbels, Ultracold 
Quantum Fields (Springer, Berlin, 2009). 

A.J. Leggett, Quantum liquids, Oxford University Press 
(New York, NY, USA, 2006). 

M. Marini, F. Pistolesi, G.C. Strinati, Eur. Phys. J. 1, 
151 (1998). 

M. Randeria, J.-M. Duan, and L.-Y. Shieh, Phys. Rev. 
B 41, 327 (1990). 

L-K. Lim, C. Morais Smith, H.T.C. Stoof, Phys. Rev. A, 
78, 013634 (2008). 

G. Bertaina and S. Giorgini, Phys. Rev. Lett. 106, 
110403 (2011). 



